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INFLUENCE OF PARENT AND DAUGHTER DECAY ON 
TRACER DIFFUSION PROFILES 
by Donald R. Behrendt 
Lewis Research Center 

SUMMARY 

In some radioactive tracer diffusion experiments, the diffusion profile is determined 
by sectioning the sample and counting the gamma activity of each of the sections. If the 
daughter decay product is radioactive and has in its decay chain a gamma ray whose * 
energy is close to that of the tracer, a correction to the counting rates of the sections 
is necessary. Because the daughter decay product is a chemical species different from 
the tracer, its diffusion coefficient will, in general, be different from the tracer. An 
equation is developed which gives the diffusion profile for the daughter for a one- 
dimensional diffusion problem. This solution can be used in making an estimate of the 
error produced by the daughter decay. With sufficient knowledge of the diffusing species, 
the correction to the diffusion profile can be made. 


INTRODUCTION 

A method commonly used to determine the diffusion profile of a radioactive tracer 
in a one dimensional diffusion experiment is to section the sample and to measure the 
activity of each of the sections by observing the decay rate of tracer to its daughter. If 
the daughter is radioactive and its activity is included with that of the parent, the result- 
ing diffusion profile will be incorrect. This problem has not been encountered very 
frequently in the past because it has been possible to choose tracers which have stable 
decay products, or, if the decay products are radioactive, to exclude by some means 
their influence on the counting rate due to the tracer. More recently, work on isotope 
effects has been done in which the dependence of diffusion on the mass of the diffusing 
atom is studied. This requires pairs of radioactive tracers of the same chemical species 
which can be diffused simultaneously into a specimen and distinguished by the charac- 
teristics of their decay. This more stringent requirement may necessitate the use of 
tracers that decay into radioactive daughters. 

An example where this occurs is the self diffusion of zirconium using zirconium 95 


as a tracer. The decay product of zirconium 95 is niobium 95 which emits a gamma ray 
whose energy is close to that of zirconium 95. 


CALCULATION 


Let n(x,t) and n^(x, t) be, respectively, the volume concentrations of the parent 
(the tracer) and the daughter decay product as a function of the distance x and time t 
where x is measured from the surface of the semi -infinite sample. One dimensional 
diffusion is considered so that the concentrations np(x, t) and n ( j(x, t) do not depend on 
y and z. Let the parent have a decay constant X p and a diffusion coefficient D p ; 
similarly, the daughter has a decay constant X^ and a diffusion coefficient D^. The 
concentration n (x, t) of the parent then satisfies the differential equation 
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ax^ 


(i) 


This is the usual one dimensional diffusion equation with the addition of the last term 
which takes into account the rate of decay of the parent. If the parent is considered to 
have an initial concentration n p Q in a narrow region of thickness a at the surface of 
the sample, the boundary conditions are 


V x ’ 0)=n pO 
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n (x, 0) = 0 a < x 
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(2b) 


and 


lim n (x, t) = 0 (2c) 

x— oo p 


From equation (2a) initially, the amount of the parent per unit area of the sample surface 



where N p has the dimensions of atoms per unit area. If the diffusion distance is suffi- 
ciently large compared with the thickness a of the plating of parent atoms, an approxi- 

* 

mation is justified in which a is allowed to approach zero in such a manner that N 
remains constant. With this approximation, the solution to equation (1) is obtained in the 
* appendix as 


n p (x,t) = 


N p exp 




(7rD p t) 1/2 


(4) 


Now, calculating the concentration of the daughter, the differential equation 
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(5) 


applies, where X p np(x,t) is the rate of production of the daughter from the decay of the 
parent, and A d n d (x,t) is the rate of loss of the daughter due to radioactive decay. 

The boundary conditions are 


n d (x, 0) = n dQ 0 < x < a 


n d (x, 0) = 0 a < x 
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and 


lim n d (x,t) = 0 
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(6b) 


(6c) 


where n d Q is the initial concentration of the daughter in the surface plating. Letting 
N d be the initial amount of the daughter per unit area of the sample surface, and, as 
before, letting a— 0 and holding N d constant where 

N d = n dO a (?) 

yields the solution to equation (5) (see appendix) as 
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The first term of equation (8) gives the contribution due to the initial surface concentra- 
tion of the daughter, and the last term gives the contribution due to the decay of the 
parent to its daughter during the diffusion anneal itself. 

The total counting rate including both parent and daughter is now calculated. The 
counting rate r p ^ of a thin section of thickness A l and area A where the center of 
the section is a distance x from the surface is 

r P , d = ” a “[ i pW x ' i > * f d x a n d (x ’ ‘>! < 9 > 

where it is assumed that the parent and daughter gamma rays cannot be distinguished. 

In equation (9), 77 represents the efficiency of counter for counting the decays and f 
and f^ are the fractions of decays of the parent and daughter, respectively, which 
result in the interfering gamma rays. 

RESULTS 

Using a digital computer and numerical integration, the integral in equation (8) was 
evaluated for various values of the and using equation (9), the total counting rate 
was evaluated. The decay constants, and A^, of both parent and daughter were set 
equal to 1x10"^ reciprocal second, a half life of about 80 days, and the diffusion time t 
was set equal to 2x10 second, about 23 days. This anneal time was chosen to be com- 
parable with the isotope half life so that the influence of the daughter on the parent dif- 
fusion profile could be clearly shown. The diffusion coefficient for the parent was 


4 



_Q 

Dp = 1x10 centimeter squared per second, and, because an examination of cases where 

the diffusion coefficient of the daughter was both larger and smaller than that of the 

-6 

parent was desired, several values of D H between 10 centimeter squared per second 
-12 u 

to 10 centimeter squared per second were used. Figure 1 which is a plot of the log of 

O 

* r p ^ ^ as a f unct * on X A> shows these results. The short-dashed curves in 
these figures represent the usual diffusion profile obtained from the counting rate of the 
parent only. Each of the other curves then represents the diffusion profile obtained by 
summing the counting rates of the parent and daughter. Figure 1(b), where D p > D^, 
shows clearly that the slopes derived from the solid curves and hence the apparent diffu- 
sion coefficients are only a few percent smaller than the true value given by the recip- 
rocal slope of the dashed line representing the parent alone. However, the influence of 
the daughter is much more marked in figure 1(b) where > Dp. In the worst case 
presented, a diffusion coefficient calculated from a straight line drawn through curves 
with D^/Dp of 10 or 100 could overestimate D p by 25 percent. An error in a diffusion 
coefficient of 5 to 10 percent is of little significance in itself. However, errors arising 
from neglect of the daughter increase as the temperature of the anneal decreases due to 
the longer annealing time required to give an adequate diffusion zone. Since the errors 
are systematic and become larger with decreasing temperature, serious error could 
result in the measured activation energy for diffusion. 

The problem produced by the activity due to daughter decay becomes most serious 
in mass dependence studies where the critical measurement is the change in the relative 
concentration of two diffusing isotopes as a function of the depth of penetration. Here, 
one is attempting to measure changes in relative concentration of only a few percent. 
Ignoring the contribution of even a small quantity of daughter activity in one of the tracer 
isotopes might easily lead to such large errors in the mass effect that the results of the 
experiment would be meaningless. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 15, 1967, 

129-03-15-01-22. 
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APPENDIX - SOLUTION OF DIFFERENTIAL EQUATION 


Parent Case 


In obtaining the solution to equations (1) and (5) we have used the method of Laplace 
transforms. Equation (1) along with the boundary conditions becomes 


ayx.t) 

3t " * 


J np(x, t) 
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- Apn p (x,t) 


n p (x,0) = n p0 [l -S a (x)] 


and 


0n (0,t) 

— P = 0 

0X 


lim n (x, t) = 0 
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where the function S a (x) is defined as 

S„(x) = 0 0 < x < a 

a 

= 1 a < x (A2) 

Transforming equation (Al) on t and letting f(x, s) be the Laplace transform of 
n p (x, t) yield 


2 

sf (x, s) - n pQ [l - S a (x)] = Dp 8 f(x ’ s) - Apf(x, s) (A3) 

0x 2 

where the boundary condition at t = 0 was used. Letting g(r, s) be the Laplace trans- 
form of f(x, s) on x, equation (A3) transforms into 


6 


sg(r, s) - _H° [l - exp(-ar)] = D [r 2 g(r, s) - rf(0, s)] - A g(r, s) 


(A4) 


. where the boundary condition at x = 0 was used. 

Equation (A4) can then be solved for g(rp d , s) and using standard methods the 
inverse Laplace transform on Tp d is obtained as 


f(x, s) = 


N p exp 

’ / s + V 

L\ D P . 

V 1/2 ' 
) 1 


_D p (s + A p )] 

1/2 


(A5) 


where a— 0 such that 


“pO'Np 


(A6) 


and Np is the surface concentration of the parent. Also, in obtaining equation (A5), 
f (0, s) has been eliminated by using the boundary condition for x-<*>. Obtaining the in- 
verse transform of f(x, s) on s in equation (A5) yields 


N P ex p - v - 


n p (x,t) = 


4D pt > 


(ffDpt) 


1/2 


(A7) 


The effect of the parent decay on the concentration is to multiply the usual result by a 
factor exp(-Apt). 


DAUGHTER CASE 


The differential equation for the concentration of the daughter is from equation (5). 


3n d (x,t) 9 2 n d (x,t) 


at 


ax' 


+ A pV X ’ *) " X d n d (x ’ 
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n d (x ’ 0) = n do[ 1 “ S a< x )] 


0n d (O,t) 

— = 0 
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and 


lim n d (x, t) = 0 

X— oo 


(A8) 


The Laplace transform on t of (A8), where f(x, s) is the transform of n d (x,t), is 


sf(x, s) - n dQ [l - S a (x)] 


/s + x\ 

4 «■ [»,C * vf'’ 


- A d f(x, s) (A9) 


where the boundary condition at t = 0 has been used. Defining g(r, s) as the Laplace 
transform of f(x, s) on s yields for equation (A9) 


sg(r, s) --^2. [l - exp(-ar)] = D d [r 2 g(r, s) - rf(0, s)] 
r 


Vp 




- * d g(r, s) (A10) 


where the boundary condition at x = 0 was used. By letting a— 0 in equation (A10) in 
such a manner that 


an do = N d (ah) 

where N d is the surface concentration of the daughter at t = 0, equation (A10) can be 
solved for g(r, s). Taking the inverse transform on r and eliminating f(0, s) by using 
the boundary condition for x— ■» yield 
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